Ab initio evaluation of the charge-ordering in a'Nal^Os 
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We report ab initio calculations of the charge ordering in a'AfctV^Os using large configurations 
interaction methods on embedded fragments. Our major result is that the 2p y electrons of the 
bridging oxygen of the rungs present a very strong magnetic character and should thus be explicitly 
considered in any relevant effective model. The most striking consequence of this result is that the 
spin and charge ordering differ substantially, as differ the experimental results depending on whether 
they are sensitive to the spin or charge density. 



The problem of charge ordering in the a'Nal^O^ com- 
pound has been widely discussed since Isobe et al |l]] 
discovered in 1996 that this compound presents a spin- 
Peierls transition at 34K. From the first X-ray structural 
determination by Carpy and Galy B a charge ordering 
was assumed in the compound. Indeed, the P1\mn non- 
centrosymmetric group yielding two inequivalent vana- 
dium atoms, it supported the hypothesis of an elec- 
tronic structure formed by an alternation of magnetic 
V i+ chains and non-magnetic V 5+ chains. When a spin- 
Peierls transition was found its interpretation through 
the dimerization of the magnetic chains seemed quite ob- 
vious. In 1998, the reinvestigation of the high temper- 
ature (HT) crystal structure || yielded a Pmmn cen- 
trosymmetric group for which the equivalence of all vana- 
dium atoms lifted any grounding to the charge ordering. 
The controversy was raised again by Damascelli et al. 
in 1998 j|. The strong on-rung doublet-doublet absorp- 
tion band, observed in the optical conductivity spectrum 
around AE — leV, was attributed by the authors to a 
charge transfer excitation of the magnetic electron be- 
tween the two vanadium atoms. This picture assumes 
a strong charge order on the rung (of about 0.8 elec- 
tron), however both the crystallographic data and the 
51 V NMR studies [[| exhibit only one type of vanadium 
atom. 

It is now widely accepted that the a'iVaV^Os spin- 
Peierls transition is in fact a double transition ||, con- 
sisting in a magneto-distortion of the lattice following 
a charge ordering of the unpaired electrons. X-ray |7]] 
and electron diffraction J^] experiments leaded to propose 
for the low temperature (LT) phase a Fmmfl symmetry 
group associated with a doubling of the unit cell in the 
a and 6 directions while it is quadrupled in the c direc- 
tion. This low temperature symmetry yields two inequiv- 
alent, but very similar, V2O5 layers, each of them con- 
taining three inequivalent vanadium atoms. 51 V NMR 
experiments ^ see however only two inequivalent sets 
of vanadium atoms. In order to explain these experi- 
mental results, a large number of hypotheses have been 



proposed for the electronic structure of the low temper- 
ature phase. Among these hypotheses the partisans of a 
total charge ordering (V 5+ and V 4+ sites) are the most 
numerous. They suppose cither the same alternation of 
magnetic V 4+ and diamagnetic V b+ chains that was first 
suggested for the high temperature phase, or a zig-zag 
charge ordering along the ladders ^| . More in agreement 
with the X-ray structural symmetries, but not with the 
NMR results, is the hypothesis of two alternated ladders, 
one supporting two equivalent l/ 4 5+ atoms and the other 
a zig-zag V 5+ -V 4+ charge ordering. Fewer authors sug- 
gest that the charge ordering is only partial [[To) and that 
one should consider the previous charge localizations hy- 
potheses with a V 4 ' 5±s , S <C 0.5 charge ordering. 

The aim of the present ab-initio quantum chemical cal- 
culations is to investigate the local electronic structure of 
the compound in order to better understand — and solve 
some of the controversies of — the charge ordering both 
in high and low temperature phases. 

For this purpose we used a method of ab-initio spec- 
troscopy of embedded crystal fragments, using large 
atomic basis sets [jllj and large differential configura- 
tions interactions methods dedicated to the accurate 
evaluation of both transitions energies and states wave- 
functions (12) (for a detailed description of the method 
applied to this system see our paper on the evaluation 
of t — J and two-band Hubbard models for the high 
temperature phase (R|]). This method has been success- 
fully used both in molecular magnetic systems jl4| and in 
strongly correlated materials p5[ | in order to determine 
super-exchange or hopping integrals with an exceptional 
accuracy. The bath is designed in order to model the 
major effects (Madelung and Pauli exclusion) of the rest 
of the crystal on the fragment, that is in such a way 
that the crystal density of states projected on the atomic 
orbitals (the most crucial being the magnetic and bridg- 
ing orbitals) is accurately reproduced in the embedded 
fragment. This goal is reached through a combination of 
punctual charges and total-ions pseudo-potentials (TIP) 
(see ref. jl3| for more details and validity checks on the 
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high temperature a'NaSfoOs phase). 

We have computed all (fundamental-)doublet to (first 
excited-) doublet excitation energies and associated wave- 
functions, on all non-equivalent pairs of nearest neigh- 
bors vanadium atoms, in the high and low temperature 
phases. The crystal fragments are composed of two adja- 
cent VO5 pyramids (corner sharing for the rungs and 
along the ladders, edge-sharing between the ladders). 
The positions of the atoms as well as the bath charges and 
TIP are taken from the X-ray crystallographic data 

The most important results of our calculations, both 
in the high and low temperature phases is that 

• the magnetic sites are not the vanadium atoms by 
themselves but the whole V — O — V rungs (as first 
suggested by Horsch [|l6| ), 

• these magnetic sites do not support one delocalized 
or partly delocalized magnetic 3d electron but that 
the 2p electrons of the bridging oxygen atoms - 
supported by the 2p y orbital pointing along the b 
axis — present also a strong magnetic character. 

Indeed the computed wave-functions (reported below for 
the HT and LT fundamental doublets) for the funda- 
mental, as well as first excited, doublet states exhibit a 
strong contribution on the configurations where the 2p y 
orbital of the bridging oxygen is half-filled. In fact the 
magnetic problem on the rungs should not be considered 
as one magnetic electron delocalized between the 3d xy of 
the two vanadium atoms but rather as a resonance be- 
tween this configuration (a\dipp) + f3\ppdj) where di and 
dj respectively stand for the magnetic 3d xy orbitals of 
the two vanadium atoms and p for the 2p y orbital of the 
bridging oxygen atom) and the configuration with three 
magnetic electrons respectively localized on the magnetic 
3d xy orbitals of two vanadiums and the 2p y orbital of the 
bridging oxygen (- (7 + S) \dipd 3 ) + ~/\d i pd j ) +8\dipdj)). 

The wave function of the HT fundamental doublet of 
the rungs is 

+0.63(|d i j>p) + |j>pri j ))/\/2 

+0.68 (-2\dipdj) + \dipdj) + \d lP dj)) /V6 

+0.13 (\dididj) + \didjdj)) /V2 

+0.12 (\di d tP ) - \pdjdj)) /V2 

+small terms . . . 

The LT phase is known for presenting an alterna- 
tion (along the a direction) between modulated and non- 
modulated ladders as well as two slightly different (a, b) 
planes. While one out of two ladders remains of simi- 
lar geometry as in the high temperature phase, the other 
ladder takes a small zig-zag form along the b direction. 
The results on the non-modulated rungs are very similar 
as in the HT phase and will not be further detailed. The 



wave functions of the LT fundamental doublets on the 
modulated rungs are 

+0.25(0. 23)\dipp) + 0.59(0.60)1^-) 
+0.53(0.52)1^^) - 0.40(0.41)1^^) - 0.13(0.11)|d»pd,-) 
+0.11(0. 11)1^^) + 0.06(0.06)1^^^) 
-0.06(0.05)1^23) + 0.09(0.09)1^4,) 
+small terms . . . 

where the number in parentheses correspond to the al- 
ternate (a,b) plane. 

The first consequence of the strong magnetic character 
of the bridging oxygen atom is that it should be explic- 
itly considered while modeling the compound. One of 
the most striking implication is that while in the single 
magnetic electron picture the charge and spin sites pop- 
ulations are equal when 3 magnetic electrons are consid- 
ered the charge and spin sites populations can strongly 
differ. Table | displays the charge and spin population of 
the magnetic orbitals. 





Charge populations 


Spin populations 




di p dj 


di p dj 


HT phase 
LT non Mod. 
LT Mod. 


0.785 1.43 0.785 
0.78 1.44 0.78 
0.61 1.45 0.94 


0.58 -0.16 0.58 
0.575 -0.15 0.575 
0.22 -0.10 0.88 



TABLE I. Magnetic orbitals charge and spin populations 
on the rungs (in electrons units) . Note that on the modulated 
rungs the Vi-O distance is shorter than the Vj—O distance. 



2 



One should notice that while the 3d charge population 
ordering on the modulated rungs is only 2S c h = 0.33, the 
participation of the oxygen to magnetic properties re- 
sults in a spin population ordering nearly twice as large 
with 2S sp = 0.66. On the non-modulated rungs how- 
ever, as expected from the crystallographic data, there is 
not any charge ordering and the overall results are very 
similar to the high temperature phase. One should how- 
ever notice that the charge and spin populations remain 
substantially different, with a strong increase of the 3d 
population (compared to the 0.5 value expected in the 
one magnetic electron picture) and a smaller increase of 
the spin polarization of the 3d orbitals by inverse polar- 
ization on the central oxygen 2p orbital. 

Another result worth to be pointed at is the substan- 
tial differences we found between the magnetic orbitals 
charge populations and the total atomic populations, ft 
is well known that such total population calculations 
are partly dependent to the choice in the atomic basis 
set, and in particular to the spatial extension of the lat- 
ter. It is however widely accepted in quantum chemistry 
that the population variations are significant and reliable 
evaluations of charge transfer processes. Table || there- 
fore presents the computed atomic populations variations 
through the transition. 





Non- modulated 


rungs 


Modulated run 


gs 




Vi O 


Vi 


Vi O 


Vi 


Charge 


-0.005 +0.000 


-0.005 


+0.007 -0.001 


-0.038 


Spin 


-0.002 +0.005 


-0.002 


-0.37 +0.05 


+0.32 



TABLE II. Variation of the total atomic charge and spin 
populations (in number of electrons units) in the low temper- 
ature phase with respect with the high temperature phase. 
Note that on the modulated rungs the Vi-O distance is shorter 



One notices immediately that while the atomic total 
spin population results are totally coherent with the mag- 
netic orbital spin populations reported in table Q, the 
atomic total charge populations are substantially differ- 
ent from the magnetic orbital charge populations. In- 
deed while the total charge disproportion between the 
two vanadium atoms of the modulated rungs is extremely 
small (2S c htot — 0.045), the disproportion between the 
two 3d magnetic orbital populations is 25 c h = 0.33. 
A finer analysis of the electrons repartition in term of 
atomic orbitals shows that there are substantial elec- 
tron population transfers between the vanadium mag- 
netic 3d xy orbitals and the other 3d vanadium orbitals. 

Indeed, while the Vj atom sees its total number of elec- 
tron decrease, the occupation of its 3d xy atomic magnetic 
orbital increases by 0.143 electron due to a population 
transfer from the 3d x 2_ y 2, 3d z 2 and 3d xz orbitals. On 
the Vi atom the population transfer goes in reverse or- 
der and the atomic magnetic orbital looses 0.166 electron 
in favor of the 3d x 2_ y 2, 3d 2 2 and 3d xz orbitals (see ta- 



ble g, |yd 



d x 



d r 



0.309 



-0.136 



-0.088 



-0.085 



TABLE III. Orbital population difference between the Vj 
and Vi atoms of the modulated rungs in the low temperature 
phase. Units are in number of electron. 



Atom 


dxy 


d>x 2 — y2 




d z 2 


Vi 


-0.166 


+0.068 


+0.044 


+0.047 


V 


+0.143 


-0.068 


-0.044 


-0.038 



TABLE IV. Difference between the 3d orbitals populations 
of the two vanadiums of the rungs in the LT phase (modulated 
rungs) and the HT phase. 



than the Vj-O distance. 
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This electron transfer can be understood if one consid- 
ers that the 2p x doubly-occupied orbital of the bridging O 
atoms strongly overlaps with the 3d x 2_ y 2, 3d z 2 and 3d xz 
orbitals of the vanadiums. V. being closer to the oxygen 
than Vj by 0.13A, the larger overlap between the oxygen 
2p x and the vanadium 3d x 2_ y 2, 3d z 2 and Zd xz orbitals in- 
creases the occupation of the later orbitals. The opposite 
phenomenum acts on the Vj atom. Such an increase in 
the 3c? atomic population would result in a proportional 
increase of the on-site repulsion unless another valence 
orbital, close to the Fermi level, is simultaneously de- 
populated. This is precisely the case with the magnetic 
3d xy orbital. The combined result of this d —> d trans- 
fers and the magnetic character of the bridging oxygen 
results in a strong difference between the total charge 
(2S c htot = 0.05) and the spin (2<5 sptot = 0.69) orderings. 

The consequence of this analysis, and this is one of 
the major conclusions of the present work to draw this 
unexpected result, is that experiments sensitive to the 
local charge density and experiments sensitive to the local 
spin density should yield substantially different results. 

Summarizing our results on the charge ordering, we 
have found (in agreement with the crystallographic data) 
that there is neither charge nor spin ordering in the high 
temperature phase, that there is neither charge nor spin 
ordering on the non-modulated rungs in the low temper- 
ature phase and that while the total charge ordering is 
very small on the modulated rungs, the spin ordering is 
on the contrary quite large. 

In the HT phase, our results agree nicely both with 
the 51 V NMR studies Q (as them we see only one 
type of vanadium atom) and optical conductivity mea- 
surements ||. Indeed, our computed on- rung doublet- 
doublet excitation energy is 1.08eV which compares very 
well with the leV optical conductivity peak. 

In the LT phase Presura et al. JTtJ conclude from op- 
tical conductivity experiments, that there should not be 
any major charge redistribution through the phase tran- 
sition, in full agreement with our results. The apparent 
contradiction between these results and the inelastic neu- 
tron scattering experiments jl^] that predict large charge 
ordering can also be lifted. Indeed, the latter are sensi- 
tive to the spin ordering which we have seen to be very 
different from the charge ordering, due to the magnetic 
nature of the bridging oxygen of the rung. Finally, the 
51 V NMR experiments || predict only two different types 
of vanadium atoms in the LT phase. NMR experiments 
are sensitive to the electronic spin polarization at the 
vanadium nuclei, that is essentially the spin polarization 
of the vanadium s orbitals induced by the spin polariza- 
tion of the magnetic sites. A simple energetic pertur- 
bative argument shows that the major contribution to 
this s spin polarization must be supported by the vana- 
dium 4s orbitals. Indeed, the 4s orbitals are the only 
ones close enough to the Fermi level to interact with the 
magnetic electrons in a non negligible way. We therefore 



computed the induced spin polarization on the vanadi- 
ums 4s orbitals. On modulated rungs we found a spin 4s 
population of 0.24 x 10~ 2 on V and 1.51 x 10~ 2 on Vj 
while on the non-modulated rungs we found 0.005 x 10~ 2 , 
that is two order of magnitude lower. As a title of com- 
parison the spin vanadium 4s polarization is in the HT 
phase 0.9 x 10~ 2 . It seems therefore possible to explain 
the NMR experiments in the low temperature phase by 
the fact that the induced spin polarization on the s or- 
bital of the non-modulated rungs vanadium atoms is so 
low that their signal is not observed in the experimental 
window. 

The main message we would like to send as a conclusion 
is that the a'AfaV^Os compound should not so much be 
considered as a quarter-filled system with one magnetic 
electron delocalized on the two vanadiums of each lad- 
der rung, but rather as a resonating system where the 
V — O — V rungs support three magnetic electrons on 
three magnetic orbitals, namely the two d x y orbitals of 
the vanadium atoms and the p y orbital of the bridging 
oxygen. A major consequence of this result is that, on 
the LT-modulated rungs, the charge and spin ordering 
are different. Indeed, while the charge ordering is very 
small, the spin ordering is quite large. The direct impli- 
cation is that experiments sensitive to the spin density 
and experiments sensitive to the charge density should 
not yield similar ordering amplitudes. We therefore think 
that a reinterpretation of experimental measurements in 
view of the present results could solve most of their ap- 
parent inconsistencies. 
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